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OPERATING  CHARACTERISTICS  FOR  INDICATOR  OR-ING 
OF  INCOHERENTLY  COMBINED  MATCHED-FILTER  OUTPUTS 

INTRODUCTION 

When  multiple  pulses  are  transmitted,  in  an  effort  to  detect  the 
presence  of  a  target,  the  multiple  echoes  should  be  optimally  processed  and 
combined  before  a  decision  is  reached.  For  received  signals  that  are 
deterministic,  except  for  independent  random  phases  between  pulses,  the 
ideal  processing  consists  of  matched  filtering,  envelope  detection,  and 
combination  according  to  a  >(n  I^  rule  [1;  chapter  VII,  (1.7)].  Since  the 
receiver  input  signal-to-noise  ratio  must  be  known  in  order  to  apply  this 
rule,  the  slightly  suboptimum  alternative  of  combining  (adding)  squared 
envelopes  is  often  adopted  [1;  ch.  VII,  (1.12)];  this  is  the  situation  to  be 
considered  here. 

In  addition,  if  the  target  has  some  movement  in  the  radial  direction, 
causing  a  Doppler  shift  of  the  echoes,  a  search  must  be  conducted  over 
frequency  at  the  receiver,  in  order  not  to  miss  the  received  signal  energy. 
For  example,  suppose  a  series  of  M  tone  bursts  at  a  common  center  frequency 
are  transmitted  and  echoed  off  a  moving  point  target.  Since  the  received 
center  frequency  will  be  unknown,  groups  of  matched  filters  will  be 
necessary,  in  order  to  cover  the  expected  range  of  frequency  shifts.  Each 
one  of  the  possible  received  center  frequencies  that  must  be  processed  is 


called  a  channel . 


In  figure  1,  a  block  diagram  of  the  processing  in  the  n-th  channel  is 
depicted.  The  M  narrowband  filters  in  the  n-th  channel  are  indicated  by 


impulse  responses  .  They  are  followed  by  detectors  which  extract 


the  squared  envelopes  of  the  filter  outputs.  These  detector  outputs  are 


then  sampled  at  times 


ui,  ■ 


which  should  correspond  to  the  times  of 


peak  signal  at  each  filter  output.  The  sampled  outputs  are  then  added,  to 
yield  channel  output  v^. 


The  block  diagram  in  figure  1  is  not  restricted  to  a  transmitted 
sequence  of  M  tone  bursts  at  a  common  center  frequency.  In  fact,  due  to  the 
general  filter  impulse  responses  and  sampling  times  allowed,  it  encompasses 
any  sequence  of  orthogonal  deterministic  signals  transmitted  at  arbitrary 
time  delays  and  frequency  offsets,  provided  they  are  known  to  the  receiver. 
The  processor  in  figure  1  also  allows  for  unknown  time  delay  to  the  target 
range  and  unknown  frequency  shift  due  to  target  movement,  by  virtue  of  the 
sampling  times  not  being  optimum,  and  the  <'ilter  impulse  responses  not  being 
matched  to  each  received  signal  component.  An  example  is  afforded  by  the 
case  where  the  filters  are  time-delayed  and/or  frequency-shifted  versions  of 
one  another. 


h 

-nm 


(D 


h(r 


T  )  exp(i2irf  T) 
nm  nm 


(1 


corresponding  to  a  time  sequence  of  frequency-stepped  pulses;  here  h  is  the 
complex  envelope  corresponding  to  impulse  response  h  [1;  pages  65-72]. 


Filters  Detectors  iampler^  Summer 

Figure  1.  Pre-Processing  for  n-th  Channel 


Figure  2.  Indicator  Or-ing  of  N  Channels 
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Another  instance  which  is  covered  by  the  processing  indicated  in  figure 
1  is  where  the  transmitted  signal  encounters  multipath  and/or  separated 
target  highlight  structure.  For  example,  a  single  transmitted  tone  burst 
might  be  received  as  four  pulses,  due  to  two  multipaths  and  two  target 
highlights.  Thus,  the  number  M  of  filters  employed  in  figure  1  is  be 
interpreted  as  the  total  number  of  received  signal  components.  Some  results 
for  the  receiver  operating  characteristics  of  this  processor  are  given  in 
[1]  and  [2], 

When  the  processing  in  the  n-th  channel  indicated  in  figure  1  is 
completed,  the  total  of  N  channels  that  must  be  considered  is  subjected  to 
the  indicator  or-ing  depicted  in  figure  2.  Namely,  the  maximum  of  the  N 


channel  outputs  is  extracted,  along  with  its  identity, 


w  =  max  ( v.| ,  Vj, 


and  compared  with  a  fixed  threshold: 


w  <  threshold:  declare  no  signal  present 


w  >  threshold:  declare  signal  present  in  channel  n 


Thus,  there  are  two  possible  outputs  from  figure  2,  the  first  being  a 


declaration  of  no  signal  present,  and  the  second  being  a  declaration  of  a 


signal  present  along  with  an  indication  of  which  channel  contains  the 


signal.  (This  latter  information  is  useful  for  identifying  the  Doppler 


shift,  for  example,  of  a  moving  target.) 


A  false  alarm  occurs  when  output  w  in  (2)  exceeds  the  threshold,  but 
there  is  no  signal  present  at  the  input.  On  the  other  hand,  a  correct 
detection  occurs  only  when  the  signal  channel  output  exceeds  the  threshold 
and  all  the  noise  channel  outputs.  That  is,  we  insist  on  accurately 
identifying  the  signal  channel,  in  order  to  achieve  a  correct  detection. 

The  performance  characteristics  of  the  processor  combination  in  figures  1 
and  2  are  of  interest,  namely  the  false  alarm  probability  and  the 
probability  of  correct  detection,  in  terms  of  M,  the  number  of  filter 
outputs  summed,  N,  the  number  of  channels  or-ed,  and  some  signal-to-noise 
ratio  measure  at  the  receiver. 

It  should  be  observed  that  the  M  received  signal  components  have  been 
presumed  to  have  undergone  no  fading.  The  only  randomness  in  the  received 
signals  are  the  independent  random  phase  shifts  between  components.  Some 
results  on  fading  signals,  including  partial  fading  between  pulses,  are 
given  in  [3];  however,  or-ing  was  not  considered  there. 

It  is  also  assumed  that  the  individual  signal  components  are  orthogonal 

with  respect  to  each  other,  perhaps  due  to  time  separation  and/or  frequency 

shift.  That  is,  at  sampling  instant  t  ,  there  is  only  one  signal 

nm 

component  contributing,  with  all  the  other  signal  components  yielding  no 
output  at  that  filter  at  that  time. 

Ihe  processor  considered  in  this  study  has  undergone  some  analysis  in 
the  past  [4];  however,  several  significant  extensions  have  been  made  here. 
First,  a  different  definition  of  detection  probability  has  been  adopted 
here,  namely  one  which  counts  as  correct  detections  only  those  events  for 
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which  the  signal  channel  output  exceeds  both  the  threshold  and  all  the  other 
undesired  noise  channel  outputs.  Second,  results  are  extended  from  a 
sinusoidal  signal  to  arbitrary  orthogonal  deterministic  signals  and  filters, 
with  arbitrary  sampling  instants;  this  allows  for  analysis  of  the  effects  of 
filter-signal  mismatch,  Doppler  offset,  time  desynchronization,  multiple 
highlights,  etc.  Third,  a  fundamentally  different  signal-to-noise  ratio 
parameter,  d,  is  used  here  to  characterize  performance,  namely,  a  measure  of 
the  total  received  signal  energy  to  noise  spectral  density  ratio,  rather 
than  the  signal-to-noise  ratio  per  pulse  (usually  assumed  identical  for  all 
pulses);  this  allows  for  arbitrary  fractional ization  of  the  total  received 
signal  energy  into  component  pulses.  Fourth,  the  detection  probability  vs. 
false  alarm  probability  curves  are  plotted  on  normal  probability  paper  with 
total  signal-to-noise  ratio,  d,  as  a  parameter;  this  straightens  out  the 
curves,  makes  them  nearly  equi-spaced  in  d,  and  affords  easy  accurate 
interpolation  in  signal-to-noise  ratio  values.  Finally,  the  current  results 
are  extended  to  much  larger  values  of  the  number,  N,  of  or-ing  channels  and 
much  smaller  false  alarm  probabilities  Pj.;  in  particular,  values  of  N  up 
to  1000,  and  values  of  Pp  as  small  as  lE-10,  are  considered. 


•A'lfmlAlk-, 
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STATISTICS  OF  FILTER  ENVELOPE-SQUARED  OUTPUT 

In  this  section,  we  derive  the  statistical  properties  of  the  output  of 
figure  1.  Suppose  a  real  narrowband  deterministic  signal  s(t)  and  a  real 
random  noise  process  n(t)  have  complex  envelopes  s(t)  and  n(t), 
respectively.  Let  the  sum  of  these  two  processes  excite  a  narrowband  filter 
h(x)  with  complex  envelope  impulse  response  h(T) .  The  complex  envelope  of 
the  filter  output  at  time  t  is  proportional  to 


c(t)  =  [s(t)  -t-  n(t)]«h(t)  =  a(t)  +  ib(t)  +  x(t)  +  iy(t),  (4) 

where 


a(t)  +  ib(t)  =  I  dT  s(T)  h(t-T) 


(5) 


is  the  deterministic  signal  output,  and 


x(t)  iy(t)  = 


dT  n(x)  h(t-T) 


(6) 


is  the  random  noise  output  process.  Then  the  filter  squared-envelope  output 
at  time  t  is 


a(t)  ib(t)  +  x(t)  +  iy{t)j  ^  = 

=  [a(t)  x(t)]^  ^  [b(t)  ^  y(t)]^ 


(7) 
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I 


More  generally,  for  M  filters,  if  signal  s^(t)  excites  filter  h^(r) , 
the  m-th  filter  squared-envelope  output  at  sample  time  t^  is 


c  (t  )|  ^  =  [a  (t  )  +  X  (t  )]^  [b^(t  )  +  y„(t„)]^  for  1  <  m  <  M. 

m'  m'l  '■  m'  m'  m'  m^  ■'  '•  m'  m  m'  m'  ■'  -  - 


Sample  times  can  be  selected  arbitrarily;  each  individual  t^^  should  be 

2  2 

chosen  to  maximize  the  size  of  the  m-th  signal  output, 


If  we  sum  these  squared-envelope  filter  output  samples,  we  have  channel 


output 


V  =  2.  )|  ^  = 

^  I  m  m  1 


The  signal  and  noise  outputs,  given  in  (5)  and  (6),  apply  for  an  arbitrary 
complex  envelope  signal  s^{t)  and  filter  h^{T)  in  the  m-th  branch  of 
the  receiver.  The  instantaneous  output  signal  squared-envelope  is 

4o)  *  OmO)  ■  K't'  *  ■  11  ^  •  <'“> 

while  the  instantaneous  output  noise  squared-envelope  is 


■  h  “ 


n(T)  h„(t-T) 


Here,  we  presume  that  a  common  broadband  noise  n(t)  excites  all  the  filters 

fh  (T»  in  the  receiver  bank.  Observe  that  if  the  m-th  signal  is  subject 
*■  m  * 

to  a  random  phase  shift,  according  to  the  factor  exp(ie  ),  this  cancels 

m 

out  of  the  envelope-squared  signal  term.  Thus,  all  the  results  here  apply 
not  only  to  a  deterministic  signal,  but  also  to  one  with  an  arbitrary  phase 
shift.  However,  no  fading  of  the  received  signal  is  allowed  in  any  of  the 
current  results. 

If  the  real  input  noise  n(t)  is  white  with  double-sided  spectral  level 
watts/Hz,  then  the  correlation  of  complex  envelope  n(t)  is  [1;  ch.  II, 
(3.11)  and  (6.22)] 


n(t)  n*(t-T)  =  4N^  6(T)  =  2Nq  i(T);  (12) 

is  the  single-sided  noise  spectral  density  level  in  watts/Hz.  By  use 
of  (6),  this  results  in  average  noise  powers  for  the  m-th  components,  as 

«^(t)  =  y^(t)  -  2N^  Jdr  IVTDp  .  (13) 

We  presume  that  all  the  filters  have  the  same  level  (energy);  thus,  we  define 
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This  is  an  important  restriction;  greater  generality  is  given  in 
[2;  appendices  B  and  C] . 

We  are  now  in  position  to  employ  the  general  results  listed  in  appendix 
A,  when  the  noise  is  Gaussian.  Namely,  define,  as  in  (A-l), 


d 


2 


7 1,  ^  = 

XIJn  1„(T)  h,(t„  -T)p 


(1 


Observe  that  the  absolute  level  of  each  filter,  h^,  cancels  out  in  this 
2 

ratio.  However,  d  does  depend  on  the  scale  of  each  signal  s^  and 
inversely  on  noise  level  . 

The  maximum  value  of  each  term  in  these  ratios  is  realized  by  choosing 


the  m-th  filter  such  that  its  impulse  response 
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where  k  is  a  complex  constant  selected  to  guarantee  the  equal  energy 
requirement  in  (14),  and  is  a  delay  inserted  for  realizability,  and  by 
choosing  sample  time  t^^  equal  to  This  is  the  matched  filter  to  the 
m-th  signal,  sampled  at  the  time  of  peak  output.  Thus,  we  have,  in  the  best 
situation , 


where  is  the  received  signal  energy  in  the  m-th  real  signal  component 
s^(t),  and  E^  is  the  total  received  signal  energy  over  all  M  paths 

2 

(branches).  Additional  interpretations  of  d  are  available  in  (A-21)  et 
seq . 

2 

This  maximum  value  of  d  in  (1?)  is  realized  only  if  the  receiving 
filters  are  the  matched  filters  (16),  and  if  the  filter  outputs  are  sampled 
at  the  correct  time  instants.  More  generally,  the  generic  value  of  d  in 
(15)  allows  for  arbitrary  signals,  filters,  and  sampling  instants,  thereby 
affording  the  possibility  of  considering  losses  due  to  mismatch  and 
desynchronization.  The  signals  can  be  time-delayed  and/or  frequency-shifted 
versions  of  each  other,  if  desired,  A  more  thorough  analysis  and  comparison 
is  presented  in  [2;  appendices  B  and  C].  The  received  signals  have 
undergone  no  fading  in  any  of  these  considerations;  thus  the  current 
analysis  applies  to  a  deterministic  signal,  except  for  random  phase. 
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Reference  to  (A-2)  and  (A-6)  now  allows  us  to  state  the  exceedance 
distribution  function  of  channel  output  v  in  (9)  as 

Prob(v  >  u)  =  1  -  P  (u)  =  0  (d,/I?/a)  for  u  ^  0, 

V  n 

where  the  Qm-function  is 

?  /  2 

0)  =  Jg  (^)  \  -2 

Parameters  d  and  r  in  (18)  are  given  by  (15)  and  (14),  respectively, 
results  pertain  to  the  signal-bearing  channel;  the  noise-only  channel 
outputs  correspond  to  setting  d  =  0. 
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FALSE  ALARM  AND  DETECTION  PROBABILITIES 


The  exceedance  distribution  function  of  the  processor  output 
for  the  n-th  channel  (see  figure  1)  is  given  by  (18)  for  signal  present  in 
that  channel.  For  those  channels  with  no  signal  present,  the  exceedance 
distribution  is 

1  -  P^°\u)  =  0^(0, T)  =  E^_.,(T^/2)  for  u  >  0,  (20) 

where  we  have  let 

T  =  )pr  /a  (21) 

for  notational  convenience,  and  defined 

En(x)  =  exp(-x)  en(x) ,  (22) 

where 

n 

en(x)  =  ^  xVk:  (23) 

k=0 


is  the  partial  exponential  [5;  6.5.11], 


false  ALARM  PROBABILITY 


Since  the  noises  in  the  N  channels  subject  to  or-ing  in  figure  2  are 
presumed  independent,  the  probability  that  a  1 1  N  outputs  do  not  exceed  a 


threshold  value  u  is 


1 


(24) 


[P 


(0) 

V 


(U)] 


N 


where  cumulative  distribution  function  was  obtained  from  (20). 

The  false  alarm  probability  is  then 


Pp  =  1  -  n  -  E„_-,(T^/2)]\ 

where  we  used  ( 21 ) . 


(25) 


DETECTION  PROBABILITY 


When  signal  is  present  in  one  channel,  we  have  several  alternative 
definitions  of  a  detection  probability.  For  example,  we  could  define  the 
probability  of  signal  detection,  as  the  probability  that  the  signal 
channel  output  exceeds  threshold  u,  disregarding  the  noise  channels 
completely;  then  directly  from  (18)  and  (21), 

PjD  -  0„(d.T),  (26) 

which  is,  of  course,  independent  of  N. 

However,  it  is  possible  that  the  noise  channels  could  also  cause  a 
threshold  crossing,  even  when  the  signal  channel  does  not.  We  can  then 
define  a  probability  of  any  detection,  as  the  probability  that  any 

channel  output  exceeds  the  threshold  u.  This  quantity  is  given  by 
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"ao  =  '  ■  V“)  - 

-  1  -  [1-  E„.,(T^/2)]"''  [1  -  0„(d,T)],  (2J) 

by  use  of  (20)  and  (18).  This  is  the  case  considered  in  [4;  see  (9)  and 
(4)  ]. 

The  problem  with  this  latter  definition  is  that,  since  we  are  interested 
in  knowing  which  channel  contains  the  signal,  the  probability  contains 
some  (rare)  events  which  indicate  the  incorrect  channel  to  contain  the 
signal.  The  best  alternative  appears  to  be  to  define  the  probability  of 
correct  detection,  as  the  probability  that  the  signal  channel  output 
exceeds  the  threshold  u  and  exceeds  all  the  noise  outputs.  In  this  case, 
the  signal  will  be  detected  and  its  channel  number  correctly  indicated. 

This  probability  is  given  by 

00 

PcD  "  i'dt  p^(t)  [P^°^t)]^~\  (28) 

u 

where  probability  density  function  p^  and  cumulative  distribution  function  P^°^ 
are  given  by  (A-4)  and  (A-9),  respectively.  Substituting  these  expressions, 
letting  x  =  {T  /a,  and  using  (21).  there  follows  the  integral  result 


'll 
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CO 


f  2  2\ 

=  I  dx  X  (J)”"''  I„_-,(dx)  expi^:^j  [l-E„_^(x^/2)] 


N-1 


(29) 


I 


From  physical  reasoning  or  mathematical  manipulations,  it  follows  that 


P  <  P  <  P  for  N  >  1 . 
CO  SO  AD 


(30) 


V 

0 

s’ 

1 


For  N  =  1,  no  or-ing,  all  three  detection  probabilities  are  equal  to 
0„(d,T). 


Also,  from  (29),  since  the  bracketed  term  is  greater  than  or  equal  to 
its  value  at  x  =  T,  we  have  the  lower  bound 


Pro  ^  q„(d,T)  for  N  >  1 


(31) 


s' 

s' 


Thus  we  have  the  tight  bounds  on  the  probability  of  correct  detection: 

-  E„_t(T^/2)]^"^  Q„(d,T)  <  Pj,p  <  Q„(d,T)  .  (32) 

To  show  how  tight  these  bounds  are,  recall  the  false  alarm  probability 
in  (25),  in  order  to  express  the  bounds  as 
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N-1 

(1  -  Pp)  ^  QH(d.T)  <  <  Q„(d.T)  . 


(33) 


For  small  false  alarm  probabilities, 

N-1 

(1  -  Pp)  ^  ii;  1  -  Pp  ^  >  1-Pp  .  (34) 


leading  to 


(1  -  Pp)  OH(d.T)  <  PcD  <  0M(d,T):  (35) 

thus  the  bounds  in  (32)  are  very  tight  for  small  false  alarm  probabilities. 
This  is  very  convenient  computationally,  since  it  means  that  we  will  not 
have  to  evaluate  the  integral  in  (29)  numerically,  but  need  only  compute  the 
simpler  quantities  0^  and  .. 

n  I 


One  special  case  of  P^^  can  be  evaluated  in  closed  form:  for  d  =  0+, 

(28)  yields 

u 

■  H’  -  ■  S '’f  ■  <36 

the  latter  relation  following  from  (25).  This  relation  agrees  with  physical 


reason i ng . 
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TIGHTNESS  OF  BOUND 

To  verify  the  accuracy  afforded  by  using  the  upper  bound  0|^(d,T), 
instead  of  the  exact  result  (29)  for  a  short  comparative  study  of  the 

two  quantities  was  conducted;  the  numerical  results  are  tabulated  in 
appendix  B.  False  alarm  probabilities  near  the  values  .1,  .01,  .001,  and 
detection  probabilities  near  the  values  .5,  .9,  .99,  .999  were  considered, 
while  M  took  on  values  1,10,  and  N  took  on  values  2,10,100,1000.  These 
ranges  of  values  encompass  most  of  the  cases  of  practical  interest;  there  is 
no  need  to  consider  smaller  values,  since  the  discrepancy  is  even 
smaller  then.  It  will  be  observed  that  for  P^  <  .1  (the  only  cases 
plotted  here),  the  differences  between  the  exact  P  and  Ou(d,T)  are 
inconsequential;  in  particular,  see  figure  B-1 . 
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GRAPHICAL  RESULTS 


In  this  section,  we  plot  the  analytical  results  for  the  false  alarm 
probability  (25)  and  the  tight  upper  bound  on  the  probability  of  correct 
detection  (33);  see  (35).  The  number  of  filter  outputs  summed,  M,  ranges 
over  the  values 


M  =  1,2,3,4,5,6.7,8,9,10  =  1(1)10, 
while  the  number  of  channel  or-ed,  N,  ranges  over  the  values 


N  =  1,  10,  100,  1000. 


The  parameter,  d,  on  the  plots  is  the  generic  signal -to-noi se  ratio  defined 
by  (15),  for  general  signals  and  filters.  The  40  combinations  corresponding 
to  (37)  and  (38)  are  plotted  on  normal  probability  paper  in  figures*  3 
through  42.  Values  of  d  small  enough  to  encompass  the  (poor  quality) 
operating  point  (Pp-^^Q)  =  (-01,. 5)  have  been  employed;  while  at  the 
high  quality  end,  values  of  d  extending  up  to  (Pp^P^p)  =  (lE-10,.999) 
have  been  used.  The  increment  in  d  is  .5  for  all  the  results  in  figures  3 
through  42. 


'All  the  figures  are  collected  together,  after  the  Summary  section. 


•,S 


■f 

J 

J 

'J 

I 
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It  will  be  observed  that  the  curves  are  approximately  equispaced  in 
parameter  d,  thereby  allowing  for  ready  accurate  interpolation  in  d,  given 
specified  Pp  and  Pj,^.  The  curves,  for  cases  in  which  N  =  1,  are 
virtually  straight  lines,  while  those  for  N  =  1000  have  developed 
significant  curvature;  nevertheless,  the  equispaced  nature  of  the  results 
readily  accommodates  interpolation  in  all  cases. 

From  these  results,  it  is  possible  to  extract  a  different  type  of 
performance  characteristic,  namely  the  required  values  of  d  to  achieve  a 
specified  quality  of  performance  in  terms  of  false  alarm  probability  and 
detection  probability.  In  figures  43  through  48,  these  results  are  plotted 
for  the  six  combinations  of 


M  =  1  ,2,4  with  P|,p  =  .5,  .9,  (39) 

while  N  varies  over  1(1)1000,  and  Pp  takes  on  the  values  IE-2,  lE-4,  IE-6, 
lR-8,  lE-10.  (Strictly,  only  the  cases  for  N  =  1,  10,  100,  1000  follow  from 
figures  3  through  42;  the  remaining  values  of  N  were  obtained  directly  from 
(25)  and  (33) . ) 

The  most  striking  feature  of  figures  43  through  48  is  their  slow 
increase  with  N,  the  number  of  channels  subjected  to  or-ing.  Certainly  the 
increase  in  required  d  values  was  anticipated,  since  or-ing  cannot  improve 
performance  capability;  however,  the  amount  of  increase  is  not  very 
significant.  Thus,  from  figure  43,  for  Pp  =  lE-10,  d  need  only  increase 
from  6.71  to  7.67  as  N  increases  from  1  (no  or-ing)  to  N  =  1000.  Greater 
increases  are  necessary  for  the  larger  Pp,  values. 
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SUMMARY 


It  will  be  easily  observed  from  the  graphical  results  in  figures  3 
through  42  that,  for  a  fixed  amount  of  or-ing  (fixed  N),  the  performance 
degrades  as  M  increases.  That  is,  for  specified  values  of  Pp  and  d,  the 
values  of  decrease  as  M  is  increased.  Alternatively,  to  maintain  a 
specified  performance  pair  values  of  d  must  be  increased  as 

M  increases.  This  is  due  to  the  fact  that  parameter  d  in  (15)  or  (17)  is  a 
total  (or  output)  signal -to-noise  ratio  measure  and  that  larger  M 
cof'responds  to  increased  fractional ization  of  the  received  signal  energy 
into  more  paths  or  branches.  Since  the  filter-output  combination  rule  is 
incoherent,  namely  adding  squared  envelopes,  this  f ract iona 1 i zation  cannot 
be  made  up  by  summation,  and  a  loss  occurs. 

On  the  other  hand,  if  we  were  to  add  more  paths  to  a  particular  system, 
then  both  M  and  d  would  increase.  Whether  this  results  in  an  improvement  or 
degradation  depends  on  the  relative  amount  of  additional  energy.  Particular 
cases  can  be  studied  quantitatively  by  referring  to  figures  3  through  42. 

In  addition,  programs  for  the  procedures  in  this  report  are  listed  in  BASIC 
in  appendix  C,  if  additional  cases  of  interest  to  the  reader  need  to  be 
invest igated . 

2 

The  maximum  value  of  d  is  given  by  (17)  as  2E.j./N^;  this  can  be 
realized  only  if  the  matched  filters  (16)  are  utilized  and  if  the  sampling 

times  are  properly  selected.  If  these  conditions  are  not  met,  the  value  of 

2 

d  given  by  (15)  must  be  employed.  In  any  event,  the  figures  are 
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parameterized  by  quantity  d,  regardless  of  what  filters  and  sampling  times 
are  used.  Thus  a  desired  value  of  d  for  a  mismatched  situation  will  require 
larger  signal  levels  for  in  (15)  than  the  values  indicated  by  the 

ideal,  (17).  In  this  manner,  the  degradation  caused  by  mismatch  and/or 
desynchronization  can  be  quantitatively  assessed. 

The  received  signal  was  assumed  to  have  undergone  no  fading  in  the 
current  analysis.  Extensions  to  fading  signals,  but  without  or-ing,  are 
available  in  [3].  This  latter  reference  presumed  a  fixed  threshold  for 
decision  variable  comparisons  (as  did  this  analysis  in  (3)  and  (18)); 
extensions  to  a  variable  threshold,  based  on  a  finite  sample  size 
noise-level  estimation  procedure,  are  currently  underway.  Results  on  this 
normalizer  in  a  fading  environment  will  be  reported  on  shortly  by  the  author 


arm 


Probability  of  Correct  Detection 


TR  8121 


E-IOE-9  E-8  C-7  E-6  £-5  £-4 


Probability  of  False  Rlarm 
Fiaure  G.  ROC  for  M= 1 .  M^IOOO 


0 


TR  81?] 


E-lOE-9  E-0  E-7  E-6  E-5  E-4  E-3  .01.02  .05 

Probability  of  False  Rlarm 
Figure  17.  ROC  for  ,  N=1CC 
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Probability  of  False  fllarm 
Figure  18.  ROC  for  M=4 ,  N= 1 000 
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Figure  26.  ROC  for  M=G ,  N=1000 
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Probability  of  False  Alarm 
Figure  28.  ROC  for  M=7 ,  N=1 
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Probability  of  False  Alarm 
Figure  34.  ROC  for  M=8 ,  N=1000 
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Figure  38.  ROC  for  ,  N-1000 
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APPENDIX  A.  Q„-FUNCTION  RELATIONSHIPS 

n 

Let  ^  and  ^  tie  independent  identically  distributed  Gaussian 

2 

random  variables,  each  with  zero  mean  and  common  variance  a  ,  and  let 

[a  1  ^  and  fb  1  ^  be  arbitrary  fixed  constants.  Define  "total" 

^  m'*  1  ‘  m'  1 

parameter 

a  m=l 

Chi-Squared  Variate 


We  are  interested  in  the  statistical  description  of  the  noncentral 
chi-squared  random  variable  of  2M  degrees  of  freedom, 


V 


[(X 


m 


+  a  ) 
m 


+  (y„ 


{A-2) 


We  will  only  list  results  here,  and  will  not  give  detailed  derivations. 


he  characteristic  function  of  v  is  [6;  page  11] 
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which  is  seen  to  depend  on  the  arbitrary  constants  |a^]  and  |b  ]  only 
2 

through  the  sum  d  in  (A-1).  The  probability  density  function  of  random 
variable  v  is  [7;  6.631  4] 


.  .  _  1  f T  /d  f  u  d  ^ 

^v^^^  ^M-1  (  a  2^2  2 


for  u  >  0 . 


The  cumulative  distribution  function  of  random  variable  v  is 


U 

(V  <  u)  =  P  (u)  =  r  dt  p  (t) 

V  sj_  » 


and  the  exceedance  distribution  function  is 


1  -  Py(u)  =  Qm(d,/ir/a)  for  u  >  0, 


where  the  Q(y|-function  is 


0^(«(.B)  ^  ^  exp^  • 

As  special  cases  of  (A-4)  and  (A-6),  for  d  =  0,  we  have  probability 


density  function 


/ 


(M-l):  (2a  ) 


I  — 2  ]  f  or  u  >  0 


and  exceedance  distribution  function 


e„_  for  „  >  0. 


where  [ 5 ;  6.5.11] 
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is  the  partial  exponential  and  where  we  define 

En(x)  =  exp(-x)  e^lx)  .  ( A-1 1 ) 

Returning  to  the  general  case  of  d  >  0  for  random  variable  v  again,  the 
cumulants  of  v  are 

X^{n)  =  (2J)  ^  for  n  >  1  .  (A-12) 

the  u-th  moments  are 

v”  =  d^/2)  for  „  >  -M.  (A-13) 

and  the  n-th  moments  are 

v"  =  (20^)"  n:  ^(-d^/2)  .  (A-14) 

Chi  Variate 


The  noncentral  chi  variate  of  2M  degrees  of  freedom  is 


z  =  V 


1/2 


fM 

.m= 


[(X  +  a  +  (y  b 

m  m  '■'m  m"^ 


(A-151 


Its  probability  density  function  is 


2  (do)  M-lvo 


exp 


2o 


2  2 


for  u  >  0,  (A-16) 


A  -3 
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and  its  exceedance  distribution  function  is 


1  -  P^Cu)  =  0M(d.  u/o)  for  u  >  0. 


(  A-1  7 


As  special  cases  of  (A-16)  and  (A-17),  for  d  =  0,  we  have  probability 
density  function 


2  u 


2M-1 


2  M 

(M-l)l  {2a 


exp 


for  u  >  0 


(A-18 


and  exceedance  distribution  function 


1  -  p<“'(u)  .  E„_,|  for  u  >  0, 


(A-19 


in  terms  of  the  functions  defined  in  (A-10)  and  (A-11). 


In  general,  for  d  >  0,  the  u-th  moment  of  random  variable  z  is 

z”  =  o”  2”^^  jF^(-v/2;M;-d^/2)  for  v  >  -  2M.  (A-20 

The  characteristic  function  and  cumulants  of  z  are  not  available  in  any 

compact  form. 


Special  Case 


If  the  constants  in  random  variable  v  in  (A-2),  and  in  random  variable  z 

in  I  A  - 1  S )  ,  satisfy 


A  4 
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=  A  cos  e^,.  =  A  sin  e^,.  (A-2i; 

where  are  arbitrary,  then  (A-l)  reduces  to 

=  M  a2/cj2,  (A-22; 

independent  of  the  particular  values  of  \ .  So  if  {e  }  were  random 

m  m 

variables  instead  of  constants,  the  statistics  of  v  and  z  in  (A-2)  and 
(A-15),  respectively,  would  be  unaffected.  This  conclusion  follows 
immediately  from  (A-3). 


In  this  latter  case  of  random  [e  ],  if  they  are  also  uniformly 

'  m 

distributed  over  2ir ,  it  is  sometimes  useful  to  define  an  individual  (common) 
s i gna I -to-noi se  ratio 


for  all  m, 


(A-23; 


hen  the  parameter  in  (A-22)  can  be  expressed  as 


More  genera  I  ) y ,  if 


d^  =  2  M  R, 


'^m  “  ^m  ®m>  *^m  '  Ap^,  sin  0^^,, 


(A-24: 


(A-25) 


where  [a^i]  are  arbitrary  constants,  then  (A-l)  reduces  to 


^  ^  ^  A^ 

2  m 

o  m-l 


( A-?b: 
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Again,  presuming  to  be  uniformly  distributed  random  variables  over 

2ir,  if  we  define  the  individual  component  signal-to-noise  ratios  as 


R 

m 


2 

h  /2 
m 


then  (A-26)  can  be  expressed  as 


d 


2 


M 


(A-27) 


(A-28) 


These  relations,  (A-24)  and  ;a-28),  afford  an  alternative  interpretation  of 

2 

the  "total"  parameter  d  in  terms  of  component  signal-to-noise  ratios. 


APPENDIX  B.  TABULATION  OF  P..  AND  O^Icl.T) 

L  U  n 

For  tne  eight  possible  combinations  of  M=1.10  with  N=2 , 1 0, 1 00, 1 000, 
vd  jes  of  tne  exact  value  of  and  the  approximation  afforded  by 

0  id.T;  are  tabulated  here.  An  explanation  of  table  B-1  ,  which  pertains 

M 

to  M-i ,  N=2 ,  fol lows : 

For  threshold  T  =  2.40,  the  false  alarm  probability  P^  =  .10912. 

Holding  these  values  fixed,  then  as  d  is  varied  from  2.2  to  5.4,  the 
detection  probabilities  vary  over  the  values  .5,  .9,  .99,  .999 
(approximately) .  This  case  is  covered  by  the  top  four  lines  in  table  B-i . 

^^hen  the  threshold  T  is  changed  to  3.25,  the  new  false  alarm  probaD'i’'ty 
is  Pp  =  .01015,  and  the  second  group  of  four  lines  in  table  6-1  pertains. 
This  procedure  is  continued  for  all  the  M,N  combinations,  while  Pp  ranges 
ovpr  the  values  .1,  .01,  .001  (approximately).  The  comparisons  for  sma’^-' 
p,  values  are  not  conducted  because  the  discrepancies  are  very  small,  as 
may  be  seen  by  inspection  of  the  tables. 

The  greatest  discrepancies  between  probabilities  P  and  Qu(d,T) 
orrijr  in  tables  B-3  and  B-4,  where  N-10.  These  particular  cases  are  plotted 
in  figure  B-1,  for  false  alarm  probabilities  in  the  .1  and  .01  regime.  For 
the  two  curves  labelled  by  A,  which  pertains  to  ,  N=10, 

P|.  .1.  show  a  very  slight  difference  between  the  two  probab  i  1  i  t  i  es  over 

ihe  range  (.5. .999).  The  label  B  actually  pertains  to  two  overlapping 
:ijr7ps  tor  M  l,  N-IQ,  P  ~  .01;  that  is,  the  plotted  values  for  P  and 


reel  Detect i on 
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Q  (d.T)  are  indistinguishable  at  this  level  of  false  alarm  probability 

M 

The  situation  for  C  and  D  is  exactly  similar,  except  that  in  these  latter 
cases,  we  have  M=10,  N=10. 


2.40 

.10912 


2.40 

2.2 

.50220 

.51005 

.10912 

3,6 

.91017 

91506 

4.6 

.98935 

.99062 

5.4 

.99890 

.99915 

3.25 

3.1 

.50353 

.50409 

.0101  5 

4.4 

.89997 

.90032 

5.5 

.99099 

.99106 

6.3 

.99920 

.99921 

3.89 

3.8 

.51653 

.51658 

.00104 

5.0 

.88951 

.88954 

6.2 

.99205 

.99206 

6.9 

.99905 

.99905 

Table  B-1  . 

Probabi 1 i ty 

Comparison 

for  M=1 ,  N=2 

d 

PCD 

OM(d.T) 

5.59 

3.5 

.49914 

.50646 

.10129 

5.1 

.89637 

.90120 

6.4 

.99083 

.99185 

7.3 

.99913 

.99931 

6.32 

4.6 

.51111 

.51166 

.01013 

6.1 

.90493 

.90526 

7.3 

.99132 

.99138 

8.1 

.99905 

.99906 

6.89 

5.3 

.49032 

.49036 

.00101 

6.8 

.90322 

.90325 

8.0 

.99162 

.99162 

8.8 

.99913 

.9991  3 

Table  B-2.  Probability  Comparison  for  M^IO,  N=-2 
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% 

< 

/  .  • 


V  ■  f  •  V 


•  W""w' I  V  '  '.'  ^  ->  -^  WT^ 


d 

^CO 

OM(d.T) 

T  =  3.70 

3.6 

,50535 

.51547 

Pp  =  .10105 

4.9 

.90037 

.90628 

6.0 

.99080 

.99193 

6.8 

.99914 

.99931 

T  =  4.29 

4.2 

,51067 

.51153 

Pp  =  .01003 

5.5 

.90499 

.90548 

6.6 

.99180 

.99189 

7.3 

.99901 

.99903 

T  =  4.79 

4.7 

.50637 

.50645 

Pp  =  .00104 

6.0 

.90378 

.90382 

7.1 

.99170 

.991 70 

7.8 

.99900 

.99900 

Table  B-5. 

Probabi 1 ity 

Comparison 

for  M=1 ,  N= 

d 

Pco 

OM(d.T) 

T  = 

6.72 

5.2 

.52240 

.53217 

II 

U- 

Q- 

.09974 

6.6 

.89913 

.90471 

7.8 

.99064 

.99169 

8.6 

.99895 

.9991  3 

T  = 

7.23 

5.8 

.51258 

.51344 

Pp  = 

.01033 

7.2 

.90082 

.90131 

8.4 

.99152 

.99160 

9.2 

.99913 

.99914 

T  = 

7.68 

6.4 

.52928 

.52936 

Pf  = 

.00102 

7.8 

.91140 

.91144 

8.9 

.99112 

.99113 

9.7 

.99910 

.99910 

Table  8-6.  Probability  Comparison  for  M-10,  N=100 


B-5 


=  4.2? 

4.2 

.51031 

.51962 

=  .10403 

5.5 

.90369 

.90884 

6.6 

.99141 

.99233 

7.3 

.99892 

.99909 

=  4.79 

4.7 

.50564 

.50645 

=  .01036 

6.0 

.90337 

.90382 

7.1 

.99162 

.991 70 

7.8 

.99899 

.99900 

=  5.25 

5.2 

.51839 

. 51846 

=  .00103 

6.5 

.90916 

.90920 

7.5 

.99008 

.99009 

8.3 

.99911 

.99911 

Table  B-7.  Probability  Comparison  for  M=1 ,  N=1000 


d 

PcD 

Om(ci.t) 

=  7.22 

5.8 

.50897 

.51785 

=  .10332 

7.2 

.89822 

.90320 

8.4 

,99098 

.99185 

9.2 

.99903 

.9991 7 

=  7.68 

6.3 

.49278 

.49354 

=  .01017 

7.8 

.91105 

.91144 

8.9 

.99105 

.99113 

9.7 

.99909 

.99910 

=  8.08 

6.8 

.50022 

.50029 

=  .00105 

8.2 

.90129 

.90133 

9.3 

.98975 

.98976 

10.1 

.99894 

.99894 

Table  B-8.  Probability  Comparison  for  M=10,  N^IOOO 


»  I  1^'  I ' 


■V'^^  II  •  f  "V^W*  »W;P  W-JP  ■■.1 


rwTTW 


Ik  cl.'l 


APPENDIX  C.  PROGRAM  LISTING 


10 
2  0 
20' 
4  0 
50 
6  0 
r  0 
20 
00 
1  00' 
1  1  0 
120 
1  2  0 
1  40 
1  5  0 
160 

1  70 
130 
1  00 
200 

2  1  0 
2  20 
2  20 
240 
250 
260 
270 
230 
230 
300 

3  1  0 
220 
2  2  0 
24  0 
2  5  O' 

6  0 
2: O' 
iiO 
2  0  0 

4  0  0 
4  1  0 
420 
4  2  0 
440 
450 
460 
470 
4  3  0 

4  00 

5  0  0 
5  10 
520 
5  20 
5  40 
550 
5  6  0 


(■’=10  I  nUNEER  OF  filter  OUTPUT;.  iUMMEL 

N=100'0  I  NUMBER  OF  OHmNNEL'2  OP-ED 

DIM  O'  100  ■  ,  Dc.‘.  1  :  10, 0:  2  ‘  '  THRE2H0LD  VHLUE  E. 

COM  Fi  '  1  00  ' ,  P'(ll  '•  100  ,  PdC'  1  00  ■ ,  Pd2‘  100  '  ,  Pd4  ■  1  00  ■  ,  Pd5  ■.  100  ' 

COM  Pd6(100>,Pd7':i00>,Pd3'lOO),Pd9'  10O',PdlO'  100' 

COM  Pd  1  1  1  00  ' ,  Pd  1  2''  100  '  ,  Pdr3<  1  00  ■'  ,  Pd  1  4  ■  1  OO  •  ,  Pd  1  5  ■  1  OO  ' 

COM  Pd  1  6  1  00  j  ,  Pd  1  7  '  1  00  ■ 

DOUBLE  M,N,I,J  I  INTEGER'2 

DhTh  2 ,  2 , 4 , 4 , 3 ,  2 , 4 , 5 ,  2 , 4 . 4 , 5 ,  2 , 4 , 4 , 5 ,  2 . 4 . 5 , 5 
doth  2 , 4 , 5 , 5 , 4 , 4 , 5 , 5 , 4 , 4 , 5 , 6 , 4 , 5 , 5 , 6 , 4 . 5 , 5 , 6 
PEhD  D'j':*:'  '  'BThRTING  '..''mLOEC  FOR  d 

0  =  0. 


U  =  U+ . 0  1 


Pt'  =  FNPf'.U,  M,  N 

IF  Pt';..l  THEN  120  I 
U  1  =  M  H  X  T  U  -  .  0  1  ,  .  0  1  > 

U  =  U+ . 0  1 

PF=FNPf TU, M, NT 

IF  Pf  lE-10  THEN  170  I 
U2  =  U 

D  e  1 1.4  =  u  2  -  u  1  '  1  0  0 . 

FOR  1=0  TO  100 
U  =  U  1  +  D  If  1  u  »  I 
IJi:  I  :'=IJ 

Pf  I  <=FNPFi'  IJ,  M,  N,' 

NEXT  I 
I  =  L  G  T  (  N  T 
Dd  =  Do'.'M,  I  T 

PRINTER  IS  PRT 
PRINT  f1,N,Do 
PRINTER  IS  CRT 
FOR  J=1  TO  17 

D  i  =  D I  j  +  ' ,  .1  -  1  1  ♦  .  5  I 

FOR  1=0  TO  100 
0  =  U'  I  '  I 


Pd  =  FNPd ' Df  , 
Pd  =  MIN'  Pd,  . 
IF  J=1  THEN 
IF  .1  =  2  then 
IF  .T  =  2  THEN 
IF  ;=4  THEN 
IF  .T  =  5  then 
IF  .1  =  6  THEN 
IF  J=7  THEN 
IF  J=3  THEN 
IF  J=3  THEN 
IF  J=10  then 
IF  .1=11  then 
IF  j=i2  Then 
IF  J=13  THEN 
IF  .1=14  THEN 
IF  J=15  THEN 
IF  J=16  THEN 
IF  .T=17  THEN 
ne;:t  I 

NEXT  J 


U  ,  M  • 

3933'2  ' 

Pdl ' I 4=Pd 
Pd2',  I  -.=Pd 
Pd  2'  I  ■=Pd 
Pd4' I 4=Pd 
Pd5'.  I  ,i=Pd 
P  d  6  *  I  4  =  P  d 
Pd7'  I  ;<=Pd 
PdS' I )=Pd 
Pd9',  I  '=Pd 
Pd  10  c:  I  >=Pd 
Pdl 1 c I '=Pd 
Pd  12'.  I  '=Pd 
Pdl 3c  I '=Pd 
Pd  1 4  c  I >=Pd 
Pdl5< I )=Pd 
Pdl6(  I >=Pd 
Pdl 7c  I  )=Pd 


UPPER  LIMIT  ON  Pf 


LOWER  LIMIT  ON  Pf 


TuTRL  DEFLECTION  PhRRMETER  d 

threshold 


C-1 


TR  8121 


FOR  1=0  TO  100 
P  t  ‘  I  '  =  F  H  I  r  i fj  1  I  P  f  •  I  J  j 
Pd  1  '  I  ;  =FM  I  nuph  I  Pd  1  I  ■  • 
Pd2CIj=FNIni-iphi  CPd2(I  >) 

P  d  3 1  I  ;!  =  F  N  I  n  M  p  h  1  ■  P  d  3  ^  I  .  '  > 

P  d  4  I  I  '  =  F  H  I  n ' .  '  p  h  1  (  ’  P  d  4  I  ■  ' 

P  d  5  '  I  J  =  F  H  I  r'l  M  p  h  1  •:  P  d  5  V  I  ■  ) 
Pdii  I  I  1  =FN  I  nvph  i  i  Pdt  (  I  ' 

P  d  7  I  I  j  =  F  H  I  n '  J  p  h  t  ‘  P  d  7  >.  I  ■  < 
PdO'  I  i  =FN  I  fi','ph  1  '  PdS'  I  '  ' 
Pd'?'  I  ‘=FNIn'.'phl  '  Pd'?'-I  '  ' 

Pd  1  0  I  =FM  I  nvph  1  '  Pd  10'  I  ■'  j 
Pdl  1  (.  I  .'=FNInvphi  I  Pdl  1  ■  1  ■  ' 
Pd  1  2  C  I  ;  =FN  I  nuph  i  i  Pd  1  2  I  >  J 
P  d  1  3  (  I  =  F  N  I  n  u  p  h  1  i  P  d  1  3  •:  I  > ) 
P  d  1  4  I  -I  =  F  H  I  ri  v  p  h  I  (  P  d  1 4  •:  I  >  > 
P  d  1  j  (  I  '  =  F  H  I  n  V  p  h  1  v!  P  d  1  5  '  I  >  > 
P  d  1  o  1  I  '  =  F  H  I  ri  M  p  h  1  c  P  d  1  c  '  I  >  > 
P  d  1  7  I  '  =  F  M  I  fi >'  p  h  I  i.  P d  1  7  '  I  >  > 
NE7;T  I 


NE7;T  I 
L  HLL  H 
END 


DEF  FMIf("phi  '.X  '  '  PM'i  jj,  , 

IF  ;-:=.5  THEN  RETURN  0, 

P=niN' X, 1 . -X> 

T=-L0G<P> 

T  =  SG(R<T  +  T) 

P=l.+T*(l.432788  +  T-*'.'.  1  89269  + T-^.  001308 


P  =  T-';  2. 515517  +  T+<:  . 
IF  X<.5  THEN  P=-P 
RETURN  P 
FNEND 


802S53  +  T+ , 0 1 0  328  ' 


DEF  FNP»  -  U,  DOUBLE  11,  N 
T=FNE'  .  5»U»U,  ri-l  ) 


FhL'OE  MLMPtI  PROBhBILIT 


P  t  =  1  .  -  1  1  .  -  T 
RETURr)  Pf 
FNEND 


DEF  FNPd'Di, 
Pd  =  FNOr(i '  M  ,  D£ 
RETURN  Pd 
FNEND 


IJ,  DOUBLE 


DETECTION  PR  0  EH  B  I  L  I  T 
UPPER  BOUND  ON  p.;d 


DEF  FNE'  ::,  DOUBLE  N 
DOUBLE  K 


'.p'  -  .  ' 

INTEGER 


T  =  '3  =  EXP' 
FOR  K  = 1 
T=T*X  K 
8  =  8  +  T 
NEXT  K 
RETURN  i 
Ft<END 


'•*  ' 
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1120 

DEF  FNCitu'- DOUBLE  M.RERL  H 

1  1  30 

Error  =  l . E-1  7 

■■ 

1140 

DOUBLE  Ml , J 

p 

1  1  50 

Q3= .  5ifR»R 

■  •  *■ 

1160 

04= . 5*B*B 

1  1  70 

Q5=EXP( 5*<Q3+Q4) ) 

S' 

1  180 

Q6=Q7=Q5 

1190 

M  1  =  M  -  1 

•  •**' 

12  00 

FOR  J=1  TO  Ml 

E 

1210 

07=07*04  J 

1  220 

06=06+07 

f 

12  30 

NEXT  J 

VL 

1240 

0  rii  =  0  5  *  0  6 

'  >  ' 

1250 

FOP  .T=l  TO  300 

V. 

1260 

05=05*03  J 

1270 

0  7  =  0  7  *  0  4  -  1  J  +  M  1  ' 

I  0  M  '  M  ,  B  ' 

'  INTEGEfti 


■ 

1  280 

06=06+07 

1  2 '90 

09=05*06 

1  300 

0  m  =  0  m  +  0  9 

13  10 

IF  0 ‘:^  ■■  =  E  f  r  0  r  *  0  fti  THEN  13  40 

L-O J 

1  320 

NEXT  J 

1  3^  3  0 

PRINT  ":-i00  TERMS  IN  FN0r.i«M, 

T 

1  340 

RETURN  MIN<  On, ,  1  . 

pr 

1  350 

FNEND 

K-.-' 

1  360 

1 

iUE  Fl  I  PLOT  PD  VS  PF  ON  HOPMmL  PPOBhBILITV  PmPF.P 

c  OM  Pf  ■;  »  > ,  Pd  1 V  *  J ,  Pd2  V  *  > ,  Pd3  >  ★  ■ ,  Pd4  <:  *  :■ ,  Pd5  ■:  > 

COM  PdS  ',  •»  j  ,  PdT  ,  PdS '  *  ,  Pd9  •,  »  ) ,  Pd  1 0 '  » 

COM  Pdll'*.',Pdl2'.*',Pdl3';*^,Pdl4<.*>,Pd'15‘*> 

COM  Pdl6<*>,Pdl7<*> • 

DIM  8$C30] , B$C32] 

DIM  X  1  ibe  1  *  <  1 :  30 ) ,  Y  1  ,ab«  1  f  U  :  30 ) 

DIM  X  ■:  o  o  r  d  1  :  3  0  ,  'i'  c  o  o  r  d  <  1  :  3  0  > 

DIM  X  g  f  t  d  '  .  1  :  3  0  > ,  Y  g  r  1  d  1 :  3  0  > 

DOUBLE  N, L  Ly, N  , Hm, I  '  INTEGERS 

I  *■ 

PI  =  "  Pr  Ob  ab  1  1  1  t  y  of  Fat  it  R  1  ar  tn  " 

B f  =  "  P r  O' b  ab  1  1  I  t  "  of  C o r  r  e  c  t  D 1 1  €  c  t  '  o n  " 


L  /  =  1  2 

REDIM  1  abt  I  f  '  1  :  L  'O.jrd'-  1  :  L  -  ,■ 

D  R  T  H  E  -  1  0  ,  E  -  9  ,  E  -  .3  ,  E  -  7  ,  E  -  E  ,  E  -  5  ,  E  -  4  ,  E  -  3  ,  .01,  .02,  .  0  o  ,  .  1 
RERD  Xlabtl-t'^*' 

DRTR  1  E  -  10  ,  1  E-9  ,  1  E-'S,  lE-T  ,  lE-E,  lE-5,  lE-4  ,  1  E  -  3  ,  .  O  1  .  .  02  ,  .  O' 


^V- 

1  560 

RERD 

C  O  0  f  d  ‘‘  ♦ 

1570 

1 

* 

1580 

L'i=  rs 

[;•? 

1  5 '90 

REDIM 

1  1  i  b  e-  1  J  *■  .  1  ;  L  ^  >  1  c  o  o  r  d 

■:  1  :  L  y  > 

1  6  0  0 

DRTR 

.  01  ,  .  02,  . 05,  .  1 ,  . 2,  .  3, 

.  4  ,  .  5  ,  .  6  ,  . 

7  ,  .  8  ,  .  9 

1  6  1  0 

DRT  R 

.95,  .'98,  .  '9  '9  ,  .  '9  '9  5  ,  .  '9  '9  S 

,  .  993 

1  620 

RERD 

1  1  i  b  €  1  J  ^ 

If.’-' 

1  6  'S-  0 

DRTR 

.  01 ,  . 02,  . 05,  .  1 ,  . 2,  .  3, 

•  4  ,  .  ,  . 

7  ,  .  3  ,  .  9 

r. 

1  6  4  0 

DRTR 

.  95  ,  .  '93,  .  99  ,  .  995,  .  998 

, .  999 

f/-' 

1  650 

RERD 

1  c  o  0  f  d  ■*  ) 

. 

1  6  6  0 

1 

di 

1670 

N  -  =  1  4 

16  8  0 

RED  I  M 

X  g  r  1  d  '  1  :  N 

r  f  k.l 

DRTR 

16-10,  1  E  -  '9  ,  1  E  -  8  .  1  E  -  7  , 

lE-6. lE-5, 

lE-4.  IE 

V.  * 

•  A  • 

1  '0  0 

RERD 

.g  r  1  '3  '  *  ' 

1  '  1  0 

1 

[R  812! 


I 

I 

A'- 


I 


k 


i^-  - 

A- 


^:- 


Lv-^ 

w 


1  :’2u 
1  ?  JO 
1  ^40 
I  ?50 
1  2&0 
1  ?70 

1  rso 
1  r  ■=*  0 
1  >:  0  0 
1  >:  1  0 


r<'..'=l3 
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